Smoking during pregnancy and its effect on placental weight: a Mendelian randomization study.

BACKGROUND: The causal relationship between maternal smoking in pregnancy and reduced offspring birth weight is well established and is likely due to impaired placental function. However, observational studies have given conflicting results on the association between smoking and placental weight. We aimed to estimate the causal effect of newly pregnant mothers quitting smoking on their placental weight at the time of delivery. METHODS: We used one-sample Mendelian randomization, drawing data from the Avon Longitudinal Study of Parents and Children (ALSPAC) (N = 690 to 804) and the Norwegian Mother, Father and Child Cohort Study (MoBa) (N = 4267 to 4606). The sample size depends on the smoking definition used for different analyses. The analysis was performed in pre-pregnancy smokers only, due to the specific role of the single-nucleotide polymorphism (SNP) rs1051730 (CHRNA5 - CHRNA3 - CHRNB4) in affecting smoking cessation but not initiation. RESULTS: Fixed effect meta-analysis showed a 182 g [95%CI: 29,335] higher placental weight for pre-pregnancy smoking mothers who continued smoking at the beginning of pregnancy, compared with those who stopped smoking. Using the number of cigarettes smoked per day in the first trimester as the exposure, the causal effect on placental weight was 11 g [95%CI: 1,21] per cigarette per day. Similarly, smoking at the end of pregnancy was causally associated with higher placental weight. Using the residuals of birth weight regressed on placental weight as the outcome, we showed evidence of lower offspring birth weight relative to the placental weight, both for continuing smoking at the start of pregnancy as well as continuing smoking throughout pregnancy (change in z-score birth weight adjusted for z-score placental weight: -0.8 [95%CI: -1.6,-0.1]). CONCLUSION: Our results suggest that continued smoking during pregnancy causes higher placental weights.

Polygenic scores for cardiovascular risk factors improve estimation of clinical outcomes in CCB treatment compared to pharmacogenetic variants alone.

Pharmacogenetic variants are associated with clinical outcomes during Calcium Channel Blocker (CCB) treatment, yet whether the effects are modified by genetically predicted clinical risk factors is unknown. We analyzed 32,000 UK Biobank participants treated with dihydropiridine CCBs (mean 5.9 years), including 23 pharmacogenetic variants, and calculated polygenic scores for systolic and diastolic blood pressures, body fat mass, and other patient characteristics. Outcomes included treatment discontinuation and heart failure. Pharmacogenetic variant rs10898815-A (NUMA1) increased discontinuation rates, highest in those with high polygenic scores for fat mass. The RYR3 variant rs877087 T-allele alone modestly increased heart failure risks versus non-carriers (HR:1.13, p = 0.02); in patients with high polygenic scores for fat mass, lean mass, and lipoprotein A, risks were substantially elevated (HR:1.55, p = 4 × 10-5). Incorporating polygenic scores for adiposity and lipoprotein A may improve risk estimates of key clinical outcomes in CCB treatment such as treatment discontinuation and heart failure, compared to pharmacogenetic variants alone.

Phenotype-based targeted treatment of SGLT2 inhibitors and GLP-1 receptor agonists in type 2 diabetes.

AIMS/HYPOTHESIS: A precision medicine approach in type 2 diabetes could enhance targeting specific glucose-lowering therapies to individual patients most likely to benefit. We aimed to use the recently developed Bayesian causal forest (BCF) method to develop and validate an individualised treatment selection algorithm for two major type 2 diabetes drug classes, sodium-glucose cotransporter 2 inhibitors (SGLT2i) and glucagon-like peptide-1 receptor agonists (GLP1-RA). METHODS: We designed a predictive algorithm using BCF to estimate individual-level conditional average treatment effects for 12-month glycaemic outcome (HbA1c) between SGLT2i and GLP1-RA, based on routine clinical features of 46,394 people with type 2 diabetes in primary care in England (Clinical Practice Research Datalink; 27,319 for model development, 19,075 for hold-out validation), with additional external validation in 2252 people with type 2 diabetes from Scotland (SCI-Diabetes [Tayside & Fife]). Differences in glycaemic outcome with GLP1-RA by sex seen in clinical data were replicated in clinical trial data (HARMONY programme: liraglutide [n=389] and albiglutide [n=1682]). As secondary outcomes, we evaluated the impacts of targeting therapy based on glycaemic response on weight change, tolerability and longer-term risk of new-onset microvascular complications, macrovascular complications and adverse kidney events. RESULTS: Model development identified marked heterogeneity in glycaemic response, with 4787 (17.5%) of the development cohort having a predicted HbA1c benefit >3 mmol/mol (>0.3%) with SGLT2i over GLP1-RA and 5551 (20.3%) having a predicted HbA1c benefit >3 mmol/mol with GLP1-RA over SGLT2i. Calibration was good in hold-back validation, and external validation in an independent Scottish dataset identified clear differences in glycaemic outcomes between those predicted to benefit from each therapy. Sex, with women markedly more responsive to GLP1-RA, was identified as a major treatment effect modifier in both the UK observational datasets and in clinical trial data: HARMONY-7 liraglutide (GLP1-RA): 4.4 mmol/mol (95% credible interval [95% CrI] 2.2, 6.3) (0.4% [95% CrI 0.2, 0.6]) greater response in women than men. Targeting the two therapies based on predicted glycaemic response was also associated with improvements in short-term tolerability and long-term risk of new-onset microvascular complications. CONCLUSIONS/INTERPRETATION: Precision medicine approaches can facilitate effective individualised treatment choice between SGLT2i and GLP1-RA therapies, and the use of routinely collected clinical features for treatment selection could support low-cost deployment in many countries.

Robust use of phenotypic heterogeneity at drug target genes for mechanistic insights: Application of cis-multivariable Mendelian randomization to GLP1R gene region.

Phenotypic heterogeneity at genomic loci encoding drug targets can be exploited by multivariable Mendelian randomization to provide insight into the pathways by which pharmacological interventions may affect disease risk. However, statistical inference in such investigations may be poor if overdispersion heterogeneity in measured genetic associations is unaccounted for. In this work, we first develop conditional F statistics for dimension-reduced genetic associations that enable more accurate measurement of phenotypic heterogeneity. We then develop a novel extension for two-sample multivariable Mendelian randomization that accounts for overdispersion heterogeneity in dimension-reduced genetic associations. Our empirical focus is to use genetic variants in the GLP1R gene region to understand the mechanism by which GLP1R agonism affects coronary artery disease (CAD) risk. Colocalization analyses indicate that distinct variants in the GLP1R gene region are associated with body mass index and type 2 diabetes (T2D). Multivariable Mendelian randomization analyses that were corrected for overdispersion heterogeneity suggest that bodyweight lowering rather than T2D liability lowering effects of GLP1R agonism are more likely contributing to reduced CAD risk. Tissue-specific analyses prioritized brain tissue as the most likely to be relevant for CAD risk, of the tissues considered. We hope the multivariable Mendelian randomization approach illustrated here is widely applicable to better understand mechanisms linking drug targets to diseases outcomes, and hence to guide drug development efforts.

Hyperglycaemia is a causal risk factor for upper limb pathologies.

BACKGROUND: Diabetes (regardless of type) and obesity are associated with a range of musculoskeletal disorders. The causal mechanisms driving these associations are unknown for many upper limb pathologies. We used genetic techniques to test the causal link between glycemia, obesity and musculoskeletal conditions. METHODS: In the UK Biobank's unrelated European cohort (N = 379 708) we performed mendelian randomisation (MR) analyses to test for a causal effect of long-term high glycaemia and adiposity on four musculoskeletal pathologies: frozen shoulder, Dupuytren's disease, carpal tunnel syndrome and trigger finger. We also performed single-gene MR using rare variants in the GCK gene. RESULTS: Using MR, we found evidence that long-term high glycaemia has a causal role in the aetiology of upper limb conditions. A 10-mmol/mol increase in genetically predicted haemoglobin A1C (HbA1c) was associated with frozen shoulder: odds ratio (OR) = 1.50 [95% confidence interval (CI), 1.20-1.88], Dupuytren's disease: OR = 1.17 (95% CI, 1.01-1.35), trigger finger: OR = 1.30 (95% CI, 1.09-1.55) and carpal tunnel syndrome: OR = 1.20 (95% CI, 1.09-1.33). Carriers of GCK mutations have increased odds of frozen shoulder: OR = 7.16 (95% CI, 2.93-17.51) and carpal tunnel syndrome: OR = 2.86 (95% CI, 1.50-5.44) but not Dupuytren's disease or trigger finger. We found evidence that an increase in genetically predicted body mass index (BMI) of 5 kg/m2 was associated with carpal tunnel syndrome: OR = 1.13 (95% CI, 1.10-1.16) and associated negatively with Dupuytren's disease: OR = 0.94 (95% CI, 0.90-0.98), but no evidence of association with frozen shoulder or trigger finger. Trigger finger (OR 1.96 (95% CI, 1.42-2.69) P = 3.6e-05) and carpal tunnel syndrome [OR 1.63 (95% CI, 1.36-1.95) P = 8.5e-08] are associated with genetically predicted unfavourable adiposity increase of one standard deviation of body fat. CONCLUSIONS: Our study consistently demonstrates a causal role of long-term high glycaemia in the aetiology of upper limb musculoskeletal conditions. Clinicians treating diabetes patients should be aware of these complications in clinic, specifically those managing the care of GCK mutation carriers. Upper limb musculoskeletal conditions should be considered diabetes complications.

Dirichlet process mixture models to impute missing predictor data in counterfactual prediction models: an application to predict optimal type 2 diabetes therapy.

BACKGROUND: The handling of missing data is a challenge for inference and regression modelling. A particular challenge is dealing with missing predictor information, particularly when trying to build and make predictions from models for use in clinical practice. METHODS: We utilise a flexible Bayesian approach for handling missing predictor information in regression models. This provides practitioners with full posterior predictive distributions for both the missing predictor information (conditional on the observed predictors) and the outcome-of-interest. We apply this approach to a previously proposed counterfactual treatment selection model for type 2 diabetes second-line therapies. Our approach combines a regression model and a Dirichlet process mixture model (DPMM), where the former defines the treatment selection model, and the latter provides a flexible way to model the joint distribution of the predictors. RESULTS: We show that DPMMs can model complex relationships between predictor variables and can provide powerful means of fitting models to incomplete data (under missing-completely-at-random and missing-at-random assumptions). This framework ensures that the posterior distribution for the parameters and the conditional average treatment effect estimates automatically reflect the additional uncertainties associated with missing data due to the hierarchical model structure. We also demonstrate that in the presence of multiple missing predictors, the DPMM model can be used to explore which variable(s), if collected, could provide the most additional information about the likely outcome. CONCLUSIONS: When developing clinical prediction models, DPMMs offer a flexible way to model complex covariate structures and handle missing predictor information. DPMM-based counterfactual prediction models can also provide additional information to support clinical decision-making, including allowing predictions with appropriate uncertainty to be made for individuals with incomplete predictor data.

Effects of physical activity and sedentary time on depression, anxiety and well-being: a bidirectional Mendelian randomisation study.

BACKGROUND: Mental health conditions represent one of the major groups of non-transmissible diseases. Physical activity (PA) and sedentary time (ST) have been shown to affect mental health outcomes in opposite directions. In this study, we use accelerometery-derived measures of PA and ST from the UK Biobank (UKB) and depression, anxiety and well-being data from the UKB mental health questionnaire as well as published summary statistics to explore the causal associations between these phenotypes. METHODS: We used MRlap to test if objectively measured PA and ST associate with mental health outcomes using UKB data and summary statistics from published genome-wide association studies. We also tested for bidirectional associations. We performed sex stratified as well as sensitivity analyses. RESULTS: Genetically instrumented higher PA was associated with lower odds of depression (OR = 0.92; 95% CI: 0.88, 0.97) and depression severity (beta = - 0.11; 95% CI: - 0.18, - 0.04), Genetically instrumented higher ST was associated higher odds of anxiety (OR = 2.59; 95% CI: 1.10, 4.60). PA was associated with higher well-being (beta = 0.11, 95% CI: 0.04; 0.18) and ST with lower well-being (beta = - 0.18; 95% CI: - 0.32, - 0.03). Similar findings were observed when stratifying by sex. There was evidence for a bidirectional relationship, with higher genetic liability to depression associated with lower PA (beta = - 0.25, 95% CI: - 0.42; - 0.08) and higher well-being associated with higher PA (beta = 0.15; 95% CI: 0.05, 0.25). CONCLUSIONS: We have demonstrated the bidirectional effects of both PA and ST on a range of mental health outcomes using objectively measured predictors and MR methods for causal inference. Our findings support a causal role for PA and ST in the development of mental health problems and in affecting well-being.

Body mass index and inflammation in depression and treatment-resistant depression: a Mendelian randomisation study.

BACKGROUND: Major depressive disorder (MDD) has a significant impact on global burden of disease. Complications in clinical management can occur when response to pharmacological modalities is considered inadequate and symptoms persist (treatment-resistant depression (TRD)). We aim to investigate inflammation, proxied by C-reactive protein (CRP) levels, and body mass index (BMI) as putative causal risk factors for depression and subsequent treatment resistance, leveraging genetic information to avoid confounding via Mendelian randomisation (MR). METHODS: We used the European UK Biobank subcohort ([Formula: see text]), the mental health questionnaire (MHQ) and clinical records. For treatment resistance, a previously curated phenotype based on general practitioner (GP) records and prescription data was employed. We applied univariable and multivariable MR models to genetically predict the exposures and assess their causal contribution to a range of depression outcomes. We used a range of univariable, multivariable and mediation MR models techniques to address our research question with maximum rigour. In addition, we developed a novel statistical procedure to apply pleiotropy-robust multivariable MR to one sample data and employed a Bayesian bootstrap procedure to accurately quantify estimate uncertainty in mediation analysis which outperforms standard approaches in sparse binary outcomes. Given the flexibility of the one-sample design, we evaluated age and sex as moderators of the effects. RESULTS: In univariable MR models, genetically predicted BMI was positively associated with depression outcomes, including MDD ([Formula: see text] ([Formula: see text] CI): 0.133(0.072, 0.205)) and TRD (0.347(0.002, 0.682)), with a larger magnitude in females and with age acting as a moderator of the effect of BMI on severity of depression (0.22(0.050, 0.389)). Multivariable MR analyses suggested an independent causal effect of BMI on TRD not through CRP (0.395(0.004, 0.732)). Our mediation analyses suggested that the effect of CRP on severity of depression was partly mediated by BMI. Individuals with TRD ([Formula: see text]) observationally had higher CRP and BMI compared with individuals with MDD alone and healthy controls. DISCUSSION: Our work supports the assertion that BMI exerts a causal effect on a range of clinical and questionnaire-based depression phenotypes, with the effect being stronger in females and in younger individuals. We show that this effect is independent of inflammation proxied by CRP levels as the effects of CRP do not persist when jointly estimated with BMI. This is consistent with previous evidence suggesting that overweight contributed to depression even in the absence of any metabolic consequences. It appears that BMI exerts an effect on TRD that persists when we account for BMI influencing MDD.

Smoking during pregnancy and its effect on placental weight: A Mendelian randomization study.

Background: The causal relationship between maternal smoking in pregnancy and reduced offspring birth weight is well established and is likely due to impaired placental function. However, observational studies have given conflicting results on the association between smoking and placental weight. We aimed to estimate the causal effect of newly pregnant mothers quitting smoking on their placental weight at the time of delivery. Methods: We used one-sample Mendelian randomization, drawing data from the Avon Longitudinal Study of Parents and Children (ALSPAC) (up to N = 805) and the Norwegian Mother, Father and Child Cohort Study (MoBa) (up to N = 4475). The analysis was performed in pre-pregnancy smokers only, due to the specific role of the genetic instrument SNP rs1051730 (CHRNA5 - CHRNA3 - CHRNB4) in affecting smoking cessation but not initiation. Results: Fixed effect meta-analysis showed a 175 g [95%CI: 16, 334] higher placental weight for pre-pregnancy smoking mothers who continued smoking at the beginning of pregnancy, compared with those who stopped smoking. Using the number of cigarettes smoked per day in the first trimester as the exposure, the causal estimate was a 12 g [95%CI: 2,22] higher placental weight per cigarette per day. Results were similar when the smoking exposures were measured at the end of pregnancy. Using the residuals of birth weight regressed on placental weight as the outcome, we showed weak evidence of lower offspring birth weight relative to the placental weight for continuing smoking. Conclusion: Our results suggest that continued smoking during pregnancy causes higher placental weights.

Association between circadian physical activity patterns and mortality in the UK Biobank.

BACKGROUND: The benefit of physical activity (PA) for increasing longevity is well-established, however, the impact of diurnal timing of PA on mortality remains poorly understood. We aimed to derive circadian PA patterns and investigate their associations with all-cause mortality. METHODS: We used 24 h PA time series from 96,351 UK Biobank participants aged between 42 and 79 years at accelerometry in 2013-2015. Functional principal component analysis (fPCA) was applied to obtain circadian PA patterns. Using multivariable Cox proportional hazard models, we related the loading scores of these fPCs to estimate risk of mortality. RESULTS: During 6.9 years of follow-up, 2,850 deaths occurred. Four distinct fPCs accounted for 96% of the variation of the accelerometry data. Using a loading score of zero (i.e., average overall PA during the day) as the reference, a fPC1 score of + 2 (high overall PA) was inversely associated with mortality (Hazard ratio, HR = 0.91; 95% CI: 0.84-0.99), whereas a score of -2 (low overall PA) was associated with higher mortality (1.69; 95% CI: 1.57-1.81; p for non-linearity < 0.001). Significant inverse linear associations with mortality were observed for engaging in midday PA instead of early and late PA (fPC3) (HR for a 1-unit increase 0.88; 95% CI: 0.83-0.93). In contrast, midday and nocturnal PA instead of early and evening PA (fPC4) were positively associated with mortality (HR for a 1-unit increase 1.16; 95% CI: 1.08-1.25). CONCLUSION: Our results suggest that it is less important during which daytime hours one is active but rather, to engage in some level of elevated PA for longevity.

BMI and well-being in people of East Asian and European ancestry: a Mendelian randomisation study.

Previous studies have linked higher body mass index (BMI) to lower subjective well-being in adult European ancestry populations. However, our understanding of these relationships across different populations is limited. Here, we investigated the association between BMI and well-being in people of (a) East Asian and (b) European ancestry in the China Kadoorie Biobank (CKB) and UK Biobank (UKB), respectively. Mendelian randomisation (MR) methods were used to test the relationship between BMI with (a) health satisfaction and (b) life satisfaction. One-sample MR enabled us to test effects in men and women separately and to test the role of cultural contexts by stratifying our analyses by urban and rural home location in both China and the UK. Further, we implemented a control function method to test the linearity of the BMI-well-being relationship. We found evidence of different associations between BMI and well-being in individuals of East Asian versus European ancestry. For example, a genetically instrumented higher BMI tentatively associated with higher health satisfaction in people of East Asian ancestry, especially in females (ß: 0.041, 95% CI: 0.002, 0.081). In contrast, there was a robust inverse association between higher genetically instrumented BMI and health satisfaction in all European ancestry UKB participants (ß: -0.183, 95% CI: -0.200, -0.165, Pdifference < 1.00E-15). We also showed the importance of considering non-linear relationships in the MR framework by providing evidence of non-linear relationships between BMI and health and life satisfaction. Overall, our study suggests potential setting-specific causality in the relationship between BMI and subjective well-being, with robust differences observed between East Asians and Europeans when considering very similar outcomes. We highlight the importance of (a) considering potential non-linear relationships in causal analyses and (b) testing causal relationships in different populations, as the casual nature of relationships, especially relationships influenced by social processes, may be setting-specific.

Mid-life leukocyte telomere length and dementia risk: An observational and mendelian randomization study of 435,046 UK Biobank participants.

Telomere attrition is one of biological aging hallmarks and may be intervened to target multiple aging-related diseases, including Alzheimer's disease and Alzheimer's disease related dementias (AD/ADRD). The objective of this study was to assess associations of leukocyte telomere length (TL) with AD/ADRD and early markers of AD/ADRD, including cognitive performance and brain magnetic resonance imaging (MRI) phenotypes. Data from European-ancestry participants in the UK Biobank (n = 435,046) were used to evaluate whether mid-life leukocyte TL is associated with incident AD/ADRD over a mean follow-up of 12.2 years. In a subsample without AD/ADRD and with brain imaging data (n = 43,390), we associated TL with brain MRI phenotypes related to AD or vascular dementia pathology. Longer TL was associated with a lower risk of incident AD/ADRD (adjusted Hazard Ratio [aHR] per SD = 0.93, 95% CI 0.90-0.96, p = 3.37 × 10-7 ). Longer TL also was associated with better cognitive performance in specific cognitive domains, larger hippocampus volume, lower total volume of white matter hyperintensities, and higher fractional anisotropy and lower mean diffusivity in the fornix. In conclusion, longer TL is inversely associated with AD/ADRD, cognitive impairment, and brain structural lesions toward the development of AD/ADRD. However, the relationships between genetically determined TL and the outcomes above were not statistically significant based on the results from Mendelian randomization analysis results. Our findings add to the literature of prioritizing risk for AD/ADRD. The causality needs to be ascertained in mechanistic studies.

Sample-constrained partial identification with application to selection bias.

Many partial identification problems can be characterized by the optimal value of a function over a set where both the function and set need to be estimated by empirical data. Despite some progress for convex problems, statistical inference in this general setting remains to be developed. To address this, we derive an asymptotically valid confidence interval for the optimal value through an appropriate relaxation of the estimated set. We then apply this general result to the problem of selection bias in population-based cohort studies. We show that existing sensitivity analyses, which are often conservative and difficult to implement, can be formulated in our framework and made significantly more informative via auxiliary information on the population. We conduct a simulation study to evaluate the finite sample performance of our inference procedure, and conclude with a substantive motivating example on the causal effect of education on income in the highly selected UK Biobank cohort. We demonstrate that our method can produce informative bounds using plausible population-level auxiliary constraints. We implement this method in the [Formula: see text] package [Formula: see text].

Sparse dimensionality reduction approaches in Mendelian randomisation with highly correlated exposures.

Multivariable Mendelian randomisation (MVMR) is an instrumental variable technique that generalises the MR framework for multiple exposures. Framed as a regression problem, it is subject to the pitfall of multicollinearity. The bias and efficiency of MVMR estimates thus depends heavily on the correlation of exposures. Dimensionality reduction techniques such as principal component analysis (PCA) provide transformations of all the included variables that are effectively uncorrelated. We propose the use of sparse PCA (sPCA) algorithms that create principal components of subsets of the exposures with the aim of providing more interpretable and reliable MR estimates. The approach consists of three steps. We first apply a sparse dimension reduction method and transform the variant-exposure summary statistics to principal components. We then choose a subset of the principal components based on data-driven cutoffs, and estimate their strength as instruments with an adjusted F-statistic. Finally, we perform MR with these transformed exposures. This pipeline is demonstrated in a simulation study of highly correlated exposures and an applied example using summary data from a genome-wide association study of 97 highly correlated lipid metabolites. As a positive control, we tested the causal associations of the transformed exposures on coronary heart disease (CHD). Compared to the conventional inverse-variance weighted MVMR method and a weak instrument robust MVMR method (MR GRAPPLE), sparse component analysis achieved a superior balance of sparsity and biologically insightful grouping of the lipid traits.

Mendelian randomisation study of body composition and depression in people of East Asian ancestry highlights potential setting-specific causality.

BACKGROUND: Extensive evidence links higher body mass index (BMI) to higher odds of depression in people of European ancestry. However, our understanding of the relationship across different settings and ancestries is limited. Here, we test the relationship between body composition and depression in people of East Asian ancestry. METHODS: Multiple Mendelian randomisation (MR) methods were used to test the relationship between (a) BMI and (b) waist-hip ratio (WHR) with depression. Firstly, we performed two-sample MR using genetic summary statistics from a recent genome-wide association study (GWAS) of depression (with 15,771 cases and 178,777 controls) in people of East Asian ancestry. We selected 838 single nucleotide polymorphisms (SNPs) correlated with BMI and 263 SNPs correlated with WHR as genetic instrumental variables to estimate the causal effect of BMI and WHR on depression using the inverse-variance weighted (IVW) method. We repeated these analyses stratifying by home location status: China versus UK or USA. Secondly, we performed one-sample MR in the China Kadoorie Biobank (CKB) in 100,377 participants. This allowed us to test the relationship separately in (a) males and females and (b) urban and rural dwellers. We also examined (c) the linearity of the BMI-depression relationship. RESULTS: Both MR analyses provided evidence that higher BMI was associated with lower odds of depression. For example, a genetically-instrumented 1-SD higher BMI in the CKB was associated with lower odds of depressive symptoms [OR: 0.77, 95% CI: 0.63, 0.95]. There was evidence of differences according to place of residence. Using the IVW method, higher BMI was associated with lower odds of depression in people of East Asian ancestry living in China but there was no evidence for an association in people of East Asian ancestry living in the USA or UK. Furthermore, higher genetic BMI was associated with differential effects in urban and rural dwellers within China. CONCLUSIONS: This study provides the first MR evidence for an inverse relationship between BMI and depression in people of East Asian ancestry. This contrasts with previous findings in European populations and therefore the public health response to obesity and depression is likely to need to differ based on sociocultural factors for example, ancestry and place of residence. This highlights the importance of setting-specific causality when using genetic causal inference approaches and data from diverse populations to test hypotheses. This is especially important when the relationship tested is not purely biological and may involve sociocultural factors.

Weak and pleiotropy robust sex-stratified Mendelian randomization in the one sample and two sample settings.

BACKGROUND: Mendelian randomization (MR) leverages genetic data as an instrumental variable to provide estimates for the causal effect of an exposure X on a health outcome Y that is robust to confounding. Unfortunately, horizontal pleiotropy-the direct association of a genetic variant with multiple phenotypes-is highly prevalent and can easily render a genetic variant an invalid instrument. METHODS: Building on existing work, we propose a simple method for leveraging sex-specific genetic associations to perform weak and pleiotropy-robust MR analysis. This is achieved by constructing an MR estimator in which pleiotropy is perfectly removed by cancellation, while placing it within the powerful machinery of the robust adjusted profile score (MR-RAPS) method. Pleiotropy cancellation has the attractive property that it removes heterogeneity and therefore justifies a statistically efficient fixed effects model. We extend the method from the typical two-sample summary-data MR setting to the one-sample setting by adapting the technique of Collider-Correction. Simulation studies and applied examples are used to assess how the sex-stratified MR-RAPS estimator performs against other common approaches. RESULTS: The sex-stratified MR-RAPS method is shown to be robust to pleiotropy even in cases where all genetic variants violated the standard Instrument Strength Independent of Direct Effect assumption. In some cases where the strength of the pleiotropic effect additionally varied by sex (and so perfect cancellation was not achieved), over-dispersed MR-RAPS implementations can still consistently estimate the true causal effect. In applied analyses, we investigate the causal effect of waist-hip ratio (WHR), an important marker of central obesity, on a range of downstream traits. While the conventional approaches suggested paradoxical links between WHR and height and body mass index, the sex-stratified approach obtained a more realistic null effect. Nonzero effects were also detected for systolic and diastolic blood pressure as well as high-density and low-density lipoprotein cholesterol. DISCUSSION: We provide a simple but attractive method for weak and pleiotropy robust causal estimation of sexually dimorphic traits on downstream outcomes, by combining several existing approaches in a novel fashion.

Investigating a possible causal relationship between maternal serum urate concentrations and offspring birthweight: a Mendelian randomization study.

BACKGROUND: Higher urate levels are associated with higher systolic blood pressure (SBP) in adults, and in pregnancy with lower offspring birthweight. Mendelian randomization (MR) analyses suggest a causal effect of higher urate on higher SBP and of higher maternal SBP on lower offspring birthweight. If urate causally reduces birthweight, it might confound the effect of SBP on birthweight. We therefore tested for a causal effect of maternal urate on offspring birthweight. METHODS: We tested the association between maternal urate levels and offspring birthweight using multivariable linear regression in the Exeter Family Study of Childhood Health (EFSOCH; n = 872) and UK Biobank (UKB; n = 133 187). We conducted two-sample MR to test for a causal effect of maternal urate [114 single-nucleotide polymorphisms (SNPs); n = 288 649 European ancestry] on offspring birthweight (n = 406 063 European ancestry; maternal SNP effect estimates adjusted for fetal effects). We assessed a causal relationship between urate and SBP using one-sample MR in UKB women (n = 199 768). RESULTS: Higher maternal urate was associated with lower offspring birthweight with similar confounder-adjusted magnitudes in EFSOCH [22 g lower birthweight per 1-SD higher urate (95% CI: -50, 6); P = 0.13] and UKB [-28 g (95% CI: -31, -25); P = 1.8 × 10-75]. The MR causal effect estimate was directionally consistent, but smaller [-11 g (95% CI: -25, 3); PIVW = 0.11]. In women, higher urate was causally associated with higher SBP [1.7 mmHg higher SBP per 1-SD higher urate (95% CI: 1.4, 2.1); P = 7.8 × 10-22], consistent with that previously published in women and men. CONCLUSION: The marked attenuation of the MR result of maternal urate on offspring birthweight compared with the multivariable regression result suggests previous observational associations may be confounded. The 95% CIs of the MR result included the null but suggest a possible small effect on birthweight. Maternal urate levels are unlikely to be an important contributor to offspring birthweight.

Calcium-channel blockers: Clinical outcome associations with reported pharmacogenetics variants in 32 000 patients.

AIMS: Pharmacogenetic variants impact dihydropyridine calcium-channel blockers (dCCBs; e.g., amlodipine) treatment efficacy, yet evidence on clinical outcomes in routine primary care is limited. Reported associations in pharmacogenomics knowledge base PharmGKB have weak supporting evidence. We aimed to estimate associations between reported pharmacogenetic variants and incident adverse events in a community-based cohort prescribed dCCB. METHODS: We analysed up to 32 360 UK Biobank participants prescribed dCCB in primary care (from UK general practices, 1990-2017). We investigated 23 genetic variants. Outcomes were incident diagnosis of coronary heart disease, heart failure (HF), chronic kidney disease, oedema and switching antihypertensive medication. RESULTS: Participants were aged 40-79 years at first dCCB prescription. Carriers of rs877087 T allele in RYR3 had increased risk of hazard ratio (HF 1.13: 95% confidence interval 1.02 to 1.25, P = .02). Although nonsignificant after multiple testing correction, the association is consistent with prior evidence. We estimated that if rs877087 T allele could experience the same treatment effect as noncarriers, the incidence of HF in patients prescribed dCCB would reduce by 9.2% (95% confidence interval 3.1 to 15.4). In patients with a history of heart disease prior to dCCB (n = 2296), rs877087 homozygotes had increased risk of new coronary heart disease or HF compared to CC variant. rs10898815 in NUMA1 and rs776746 in CYP3A5 increased likelihood of switching to an alternative antihypertensive. The remaining variants were not strongly or consistently associated with studied outcomes. CONCLUSION: Patients with common genetic variants in NUMA1, CYP3A5 and RYR3 had increased adverse clinical outcomes. Work is needed to establish whether outcomes of dCCB prescribing could be improved by prior knowledge of pharmacogenetics variants supported by clinical evidence of association with adverse events.

Islet Autoantibody Level Distribution in Type 1 Diabetes and Their Association With Genetic and Clinical Characteristics.

CONTEXT: The importance of the autoantibody level at diagnosis of type 1 diabetes (T1D) is not clear. OBJECTIVE: We aimed to assess the association of glutamate decarboxylase (GADA), islet antigen-2 (IA-2A), and zinc transporter 8 (ZnT8A) autoantibody levels with clinical and genetic characteristics at diagnosis of T1D. METHODS: We conducted a prospective, cross-sectional study. GADA, IA-2A, and ZnT8A were measured in 1644 individuals with T1D at diagnosis using radiobinding assays. Associations between autoantibody levels and the clinical and genetic characteristics for individuals were assessed in those positive for these autoantibodies. We performed replication in an independent cohort of 449 people with T1D. RESULTS: GADA and IA-2A levels exhibited a bimodal distribution at diagnosis. High GADA level was associated with older age at diagnosis (median 27 years vs 19 years, P = 9 × 10-17), female sex (52% vs 37%, P = 1 × 10-8), other autoimmune diseases (13% vs 6%, P = 3 × 10-6), and HLA-DR3-DQ2 (58% vs 51%, P = .006). High IA-2A level was associated with younger age of diagnosis (median 17 years vs 23 years, P = 3 × 10-7), HLA-DR4-DQ8 (66% vs 50%, P = 1 × 10-6), and ZnT8A positivity (77% vs 52%, P = 1 × 10-15). We replicated our findings in an independent cohort of 449 people with T1D where autoantibodies were measured using enzyme-linked immunosorbent assays. CONCLUSION: Islet autoantibody levels provide additional information over positivity in T1D at diagnosis. Bimodality of GADA and IA-2A autoantibody levels highlights the novel aspect of heterogeneity of T1D. This may have implications for T1D prediction, treatment, and pathogenesis.

Should RECOVERY have used response adaptive randomisation? Evidence from a simulation study.

BACKGROUND: The Randomised Evaluation of COVID-19 Therapy (RECOVERY) trial is aimed at addressing the urgent need to find effective treatments for patients hospitalised with suspected or confirmed COVID-19. The trial has had many successes, including discovering that dexamethasone is effective at reducing COVID-19 mortality, the first treatment to reach this milestone in a randomised controlled trial. Despite this, it continues to use standard or 'fixed' randomisation to allocate patients to treatments. We assessed the impact of implementing response adaptive randomisation within RECOVERY using an array of performance measures, to learn if it could be beneficial going forward. This design feature has recently been implemented within the REMAP-CAP platform trial. METHODS: Trial data was simulated to closely match the data for patients allocated to standard care, dexamethasone, hydroxychloroquine, or lopinavir-ritonavir in the RECOVERY trial from March-June 2020, representing four out of five arms tested throughout this period. Trials were simulated in both a two-arm trial setting using standard care and dexamethasone, and a four-arm trial setting utilising all above treatments. Two forms of fixed randomisation and two forms of response-adaptive randomisation were tested. In the two-arm setting, response-adaptive randomisation was implemented across both trial arms, whereas in the four-arm setting it was implemented in the three non-standard care arms only. In the two-arm trial, randomisation strategies were performed at the whole trial level as well as within three pre-specified patient subgroups defined by patients' respiratory support level. RESULTS: All response-adaptive randomisation strategies led to more patients being given dexamethasone and a lower mortality rate in the trial. Subgroup specific response-adaptive randomisation reduced mortality rates even further. In the two-arm trial, response-adaptive randomisation reduced statistical power compared to FR, with subgroup level adaptive randomisation exhibiting the largest power reduction. In the four-arm trial, response-adaptive randomisation increased statistical power in the dexamethasone arm but reduced statistical power in the lopinavir arm. Response-adaptive randomisation did not induce any meaningful bias in treatment effect estimates nor did it cause any inflation in the type 1 error rate. CONCLUSIONS: Using response-adaptive randomisation within RECOVERY could have increased the number of patients receiving the optimal COVID-19 treatment during the trial, while reducing the number of patients needed to attain the same study power as the original study. This would likely have reduced patient deaths during the trial and lead to dexamethasone being declared effective sooner. Deciding how to balance the needs of patients within a trial and future patients who have yet to fall ill is an important ethical question for the trials community to address. Response-adaptive randomisation deserves to be considered as a design feature in future trials of COVID-19 and other diseases.